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1.  Introduction 


Traditionally  in  aeroballistic  applications,  the  pitch-damping  coefficient  has  been  treated  as  the 
sum  of  two  individual  coefficients  that  produce  an  aerodynamic  moment  proportional  to  the 
angular  rate  associated  with  the  angle  of  attack.  In  fact,  these  two  individual  coefficients 
represent  moments  proportional  to  two  different  angular  rates,  although  for  many 
nonmaneuvering  flight  trajectories,  including  those  flown  in  ballistic  aerodynamic  ranges,  these 
two  angular  rates  are  essentially  equivalent.  For  this  reason,  the  pitch-damping  coefficient  sum 
is  often  treated  as  a  single  parameter.  In  some  cases,  such  as  for  maneuvering  flight  vehicles,  the 
simplification  that  the  two  angular  rates  are  equivalent  is  no  longer  valid,  and  the  two  individual 
components  of  the  pitch-damping  coefficient  sum  must  be  determined  independently. 
Experimental  determination  of  these  coefficients  is  difficult  at  best  and  is  probably  not  possible 
using  traditional  aerodynamic  ranges. 

Prior  research  has  focused  on  the  development  and  application  of  computational  methods  for 
predicting  the  pitch-damping  coefficient  sum.  The  basis  of  the  technique  is  to  impose  a 
particular  motion  (coning  motion)  on  the  flight  vehicle,  which  produces  moments  that  are 
proportional  to  the  pitch-damping  coefficient  sum.  Because  coning  motion  produces  a  steady 
flow  field  in  the  cases  of  interest,  efficient  numerical  techniques  can  be  readily  applied  and  the 
pitch-damping  coefficient  sum  determined  in  a  cost-efficient  manner.  These  techniques  have 
been  benchmarked  with  experimental  data,  and  good  agreement  between  computation  and 
experiment  has  been  found  for  a  number  of  flight  vehicle  geometries  {1-10). 

In  the  current  effort,  the  method  is  extended  so  that  the  individual  components  of  the  pitch¬ 
damping  coefficient  can  be  determined  independently.  The  method  uses  imposed  motions  that 
excite  the  two  angular  rates  independently  so  that  the  forces  and  moments  attributable  to  these 
rates  can  be  assessed  separately.  A  key  feature  of  the  approach  is  that  the  motions  produce 
steady  flow  fields  that  can  be  computed  in  a  cost-efficient  manner.  The  techniques  described 
here  represent  a  unique  aerodynamic  capability  for  a  problem  that  is  difficult  to  address 
experimentally.  Other  than  approximate  methods,  these  results  most  likely  represent  the  first 
numerical  predictions  of  the  individual  components  of  the  pitch-damping  coefficient  sum  for 
flight  vehicles.  Subsequent  to  its  original  publication  {8-10),  the  technique  has  also  been 
adopted  by  other  researchers  {11). 

This  report  presents  a  derivation  of  the  transverse  aerodynamic  force  and  moment  equations 
associated  with  the  motions  of  interest  using  a  general  force  and  moment  expansion  for 
symmetric  flight  bodies.  The  resulting  equations  demonstrate  that  the  aerodynamic  forces  and 
moments  can  be  excited  independently  using  the  appropriate  motions.  Results  are  presented  for 
a  family  of  axisymmetric  flight  bodies  at  supersonic  flight  velocities.  A  schematic  of  the  Army- 
Navy  spinner  rocket  (ANSR)  flight  body  geometries  is  shown  in  figure  1 .  Both  of  the 
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components  and  )  of  the  piteh-damping  coeffieient  sum,  ,  are  determined 

independently.  Comparison  of  the  predieted  aerodynamic  coefficients  is  made  with  slender  body 
theory  and  numerieal  solutions  of  the  nonlinear  unsteady  potential  equation.  The  sum  of  the  two 
individually  determined  coeffieients  is  also  compared  with  the  piteh-damping  coeffieient  sum, 

’  determined  from  eoning  motion.  Both  methods  of  determining  the  piteh-damping 

coeffieient  sum  are  in  good  agreement  with  experimental  data. 


2.  Theoretical  Background 


2,1  The  Force  and  Moment  Expansion 

ft  is  common  in  many  aerodynamics  applications  to  use  a  body-fixed,  nonrolling  coordinate 
system  to  describe  both  the  dynamics  and  the  system  of  forces  and  moments  that  act  on  the  flight 
vehicle  (12).  The  nonrolling  coordinate  system  allows  the  description  of  the  vehicle  dynamics  to 
be  simplified  for  certain  classes  of  flight  vehicles  that  possess  particular  types  of  geometric 
symmetry.  Rotationally  symmetric  flight  vehicles,  which  are  the  focus  of  the  current  research, 
represent  one  class  of  vehicles  where  the  nonrolling  frame  has  been  effectively  (and 
traditionally)  used.  For  more  complicated  geometries,  such  as  aircraft,  the  advantages  of  the 
nonrolling  frame  are  reduced  and  other  coordinate  frames  such  as  a  completely  body-fixed 
coordinate  system  are  typically  used. 
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In  the  current  effort,  the  primary  reason  for  initially  describing  the  aerodynamic  forces  and 
moments  using  the  nonrolling  coordinate  system  is  the  fact  that  the  description  is  well 
established  for  symmetric  flight  vehicles.  The  nonrolling  coordinate  frame  is  an  orthogonal 
right-handed  system  (x,  y,  z)  centered  at  the  body’s  center  of  gravity  (CG).  The  x  -axis  is 
aligned  along  the  projectile  longitudinal  axis  with  the  positive  direction  oriented  toward  the 
projectile  nose.  The  z  -axis  is  “initially”  oriented  downward  with  the  x  -  z  plane  perpendicular 
to  the  ground.  The  angular  motion  of  the  nonrolling  coordinate  frame  is  such  that,  with  respect 
to  an  inertial  frame,  the  x  -component  of  the  coordinate  frame’s  angular  velocity  is  zero. 
Although  the  time-dependent  orientation  of  the  nonrolling  frame  may  be  difficult  to  visualize, 
the  nonrolling  frame  is  essentially  equivalent  to  the  “fixed-plane”  coordinate  system  for  small 
amplitude  motions.  In  the  fixed-plane  coordinate  system,  the  x  -  z  plane  remains  perpendicular 
to  the  ground  for  all  time.  The  total  angular  velocity  of  the  flight  vehicle  can  be  described  in 
terms  of  its  angular  velocity  components  (p,  q,  7)  along  the  x  ,  y  ,  and  z  axes,  respectively. 

The  angular  velocity  of  the  nonrolling  frame  can  be  described  in  terms  of  the  transverse  angular 
velocities  q  and  v  because  the  angular  velocity  of  the  nonrolling  frame  along  the  x  -axis  is 
always  zero.  The  flight  body  may,  however,  have  a  nonzero  spin  rate,  p,  about  its  longitudinal 
axis.  Further  details  about  these  coordinate  frames  are  discussed  in  reference  (12). 


The  moment  expansion  for  a  rotationally  symmetric  missile  in  the  nonrolling  coordinate  frame  is 
shown  in  equation  1 .  This  moment  expansion  is  similar  to  the  moment  proposed  by  Murphy 
(12).  The  moment  formulation  uses  complex  variables  to  separate  the  moment  components,  Cj„ 
and  C„ ,  that  are  oriented  along  the  y  and  z  axes,  respectively.  The  third  moment  component, 
the  roll  moment,  can  be  handled  separately  and  is  not  of  consequence  in  this  study. 


C„.+iC„  = 


Cn 

-iC^ 

r'pa 

(1) 


In  the  moment  expansion,  the  pitching  moment  coefficient  slope,  ,  and  the  coefficient 

represent  moments  that  are  proportional  to  the  complex  yaw,  ^  ,  and  yawing  rate,  ^ ' , 
respectively. 


The  complex  yaw  and  yawing  rate  are  defined  as  follows.  (In  the  analysis  presented  here,  there  is 
no  need  to  distinguish  between  pitch  and  yaw,  and  the  terms  may  be  interchanged.)  The  usage 
follows  that  of  Murphy  (12). 


V  +1W 

V 


(2) 


and 


(3) 
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Here,  v  and  w  are  the  y  -  and  z  -components  of  the  velocity  vector,  ,  that  describes  the 
velocity  of  the  body  CG  relative  to  the  inertial  frame.  The  magnitude  of  this  vector  is  denoted  as 
V.  The  angular  rate  is  obtained  by  taking  the  derivative  of  the  complex  yaw  with  respect  to  the 
flight  path  coordinate,  s,  which  is  nondimensionalized  by  some  characteristic  length,  I ,  typically 
the  body  diameter. 


The  coefficient,  ,  represents  a  moment  that  is  proportional  to  the  complex  transverse  angular 
velocity  of  the  vehicle,  Jl ,  as  defined  in  the  following; 


.....  (q  +  iTV 

//  =  ^ ^ 

V 


(4) 


Here,  q  and  ?  are  the  y  -  and  z  -components  of  angular  velocity  of  the  vehicle  in  the 
nonrolling  coordinate  system.  The  remaining  coefficient  in  the  moment  expansion,  the  Magnus 
moment  coefficient,  ,  accounts  for  a  side  moment  due  to  flow  asymmetries  produced  by  the 

combination  of  spin  and  yaw. 


It  should  be  noted  that  the  moment  formulation  neglects  the  variation  of  the  moments  with  roll 
angle  under  the  assumption  that  these  variations  are  small.  For  axisymmetric  vehicles,  the 
variations  with  roll  angle  should  not  exist  because  the  geometry  will  not  change  as  the  roll 
orientation  changes.  Roll  variations  in  the  aerodynamic  coefficients  for  other  types  of 
rotationally  symmetric  vehicles  are  typically  negligible  for  small  amplitude  motions.  In  general, 
roll  variations  may  be  difficult  to  detect  in  flight  because  the  effect  of  roll  orientation  tends  to  be 
averaged  out  over  the  course  of  a  yaw  cycle  if  the  body  is  spinning. 


2,2  Planar  Motions 

For  planar  motions,  the  aerodynamic  moments  that  act  on  an  axisymmetric  flight  body  can  be 
written  in  terms  of  the  following  expansion.  The  force  expansion  has  a  similar  form. 


C„+iC  = 


i 

p 

-h 

n 

1 

sin  «  + 

( 

^  di^ 

LILq, 

niq 

(yj 

(5) 


Two  damping  moments,  and  ,  represent  moments  proportional  to  the  angular  rates  d 

and  q  that  are  associated  with  the  angle  of  attack,  a,  and  the  angular  displacement  of  the 
longitudinal  axis  of  the  body  with  respect  to  the  earth-fixed  axis  system  denoted  by  the  angle,  0 , 
respectively.  These  angles  are  displayed  schematically  in  figure  2  for  the  case  of  planar  motion. 

For  a  typical  ballistic  trajectory,  over  the  course  of  many  yaw  cycles,  the  flight  body  will  travel  a 
nearly  rectilinear  flight  path.  In  this  case,  the  free-stream  velocity  vector,  ,  has  a  fixed 
orientation  with  respect  to  the  earth-fixed  axes,  Xe  and  Ze.  For  a  rectilinear  flight  path,  the 
angular  rates  d  and  q  will  be  equal.  The  moment  expansion  can  be  simplified  by  combining  the 
two  damping  coefficients  into  a  single  coefficient  sum,  which  is  proportional  to  a  single  angular 
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rate  (either  cc  or  q).  For  maneuvering  flight  bodies,  this  simplification  may  not  be  valid,  and  the 
individual  coefficients  must  be  treated  independently. 


To  produce  aerodynamic  moments  proportional  to  the  damping  terms,  one  or  both  of  the  angular 
velocities  must  be  nonzero.  One  simple  planar  motion  that  produces  a  nonzero  q  and  zero  d 
angular  velocity  is  a  circular  looping  motion  as  shown  in  figure  3.  If  such  a  motion  is  performed 
at  constant  angular  velocity,  Q ,  (q  =  ^  =  Q)  ,  and  constant  angle  of  attack  (which  may  or  may 
not  be  zero),  the  following  form  of  the  force  and  moment  expansions  results: 


C,.  +iC„  = 


vVy 


C„  +C^ 

pa  ^^^a 


sin«  +  C, 


(6) 


Here,  the  in-plane  moment,  Cj^ ,  contains  contributions  from  both  the  pitching  moment  slope, 

,  and  the  pitch-damping  coefficient,  .  The  side  moment,  Cn,  is  identical  in  form  to  the 

side  moment  due  to  constant  a  motion  and  is  independent  of  the  angular  rate  q.  Because  the 
damping  terms  are  independent  of  the  spin  rate  and  the  side  moment  is  independent  of  the 
angular  rate  q,  the  spin  rate  is  assumed  to  be  zero  to  simplify  the  discussion  of  the  looping 
motion  that  follows. 


With  respect  to  the  inertial  frame  of  reference,  the  flow  field  is  periodic  and  unsteady.  However, 
in  the  nonrolling  frame  (shown  in  figure  3  as  the  x  -  z  axes),  the  flow  field  is  potentially  steady. 
Indeed,  the  moment  expansion  displays  no  unsteadiness  since  all  the  terms  on  the  right-hand  side 
are  constants.  Because  the  flow  field  is  potentially  steady  in  the  nonrolling  frame  for  this  type  of 
motion,  this  frame  is  suitable  for  use  with  computational  approaches  that  are  based  on  steady 
flow  techniques. 
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Figure  3.  Schematic  of  planar  looping  motion. 


For  looping  motion,  the  moment  expansion  contains  contributions  from  two  of  the  aerodynamic 
coefficients.  The  contribution  from  the  pitching  moment  coefficient  slope,  ,  can  be 

eliminated  if  the  body  longitudinal  axis  is  tangent  to  the  radius  of  curvature  of  the  loop.  In  this 
case,  the  angle  of  attack  is  zero,  and  the  damping  moment,  ,  is  directly  proportional  to  the 

net  aerodynamic  moment.  If  the  damping  moment  exhibits  a  dependence  on  angle  of  attack, 
several  experiments  or  computations  involving  different  looping  rates  must  be  performed  for 
each  angle  of  attack.  (If  the  looping  rate  is  varied,  the  loop  radius  must  also  be  changed  to 
maintain  a  constant  vehicle  velocity:  V  =  •)  The  damping  moment  can  be  obtained  by 

computing  the  variation  in  the  net  aerodynamic  moment  with  looping  rate.  If  the  net 
aerodynamic  moment  displays  a  nonlinear  variation  with  both  angle  of  attack  and  looping  rate,  it 
may  be  difficult  to  separate  the  contributions  attributable  to  the  and  ,  without 

additional  information  or  assumptions.  In  general,  variations  of  the  coefficients  with  angle  of 
attack  might  be  reasonably  expected,  while  variations  in  the  coefficients  with  angular  rate  are 
less  common.  In  any  event,  the  looping  motion  can  be  used  to  check  for  these  types  of  nonlinear 
variations. 

One  feature  of  this  type  of  motion  is  that  high  angular  velocities  or  large  loop  radii  are  required 
to  generate  high  velocities,  making  this  type  of  motion  impractical  for  use  in  experimental 
testing.  However,  this  is  not  a  problem  for  computational  approaches  because  large  loop  radii  or 
high  angular  velocities  can  be  easily  accommodated. 
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2,3  Helical  Motions 

Other  more  complicated  motions  can  be  proposed,  which  produce  aerodynamic  moments 
proportional  to  the  angular  rates,  fl  and  ^ ' ,  while  still  producing  steady  flow  fields  when 
viewed  from  the  appropriate  coordinate  frame.  Two  such  motions  require  the  CG  of  the  flight 
vehicle  to  traverse  a  helical  flight  path.  The  first  motion  requires  the  vehicle’s  longitudinal  axis 
to  be  oriented  in  the  same  direction  as  the  center  of  rotation  of  the  helix  but  displaced  by  a 
constant  distance.  Figure  4  shows  a  three-dimensional  (3-D)  view  of  the  motion.  A  two- 
dimensional  (2-D)  projection  of  this  motion  on  the  vertical  plane  is  shown  in  figure  5.  This 
particular  motion  produces  no  rotation  of  the  nonrolling  coordinate  frame  relative  to  an  earth- 
fixed  coordinate  frame,  and  hence,  the  angular  velocity  p.  is  zero.  The  angle  of  attack  and  its 

angular  rate  vary  continuously,  producing  moment  components  associated  with  the  coefficients 
Cj^  and  ,  respectively.  This  motion  is  referred  to  as  “q  =  0  helical  motion”  because  the 

angular  rates  associated  with  the  damping  coefficient  are  zero. 


For  the  second  motion,  the  longitudinal  axis  of  the  flight  vehicle  remains  tangent  to  the  helical 
flight  path  at  each  point  along  the  trajectory.  Figure  6  shows  a  3-D  view  of  this  motion.  A  2-D 
projection  of  this  motion  on  the  vertical  plane  is  also  shown  in  figure  7.  The  angle  of  attack  of 
the  incident  airstream  is  zero  because  both  the  longitudinal  axis  of  the  body  and  the  free-stream 
velocity  vector  are  tangent  to  the  flight  path.  The  resulting  yawing  rate  is  also  zero  because  the 
angle  of  attack  is  constant.  The  angular  orientation  of  the  flight  body  changes  continuously  with 
respect  to  an  earth-fixed  reference  frame,  producing  a  nonzero  angular  rate,  .  As  a  result, 
moment  components  associated  with  the  damping  moment  are  produced.  This  motion  is 

referred  to  as  “  dr  =  0  helical  motion”  because  the  angular  rates  associated  with  the  damping 
coefficient  are  zero.  The  dr  =  0  helical  motion  produces  a  steady  flow  field  when  viewed 

from  the  appropriate  coordinate  system,  while  its  2-D  projection  is  clearly  time  dependent. 

The  net  transverse  aerodynamic  moment  in  the  nonrolling  frame  can  be  determined  for  each  of 
the  motions  by  substituting  the  appropriately  determined  angle  of  attack,  ^  ,  and  angular  rates, 

^ '  and  p ,  into  equation  1 .  The  angle  of  attack  and  angular  rates  for  both  types  of  helical 
motion  are  shown  in  equations  7-14. 


q  =  0  Helical  Motion 


~  V  +  iw  „ 

C  = - =  -0  Q 

V 


int 


R„Q 


a/u'+^^'Ro 


Rq^ 

V 


(V) 

(8) 
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4=4'  =  o. 


(14) 


Here,  Q  is  the  angular  velocity  of  the  body  about  the  helix  axis,  Rq  is  the  perpendicular  distance 
between  the  helix  axis  and  the  body  CG,  U  is  the  component  of  velocity  along  the  helix  axis,  and 
V  is  the  total  linear  velocity  of  the  CG. 

To  completely  define  the  motions,  the  spin  rate  of  the  body  p  needs  to  be  specified.  It  is 
convenient  to  set  the  spin  rate  of  the  body  p  to  zero  because  this  will  eliminate  the  contribution  to 
the  net  aerodynamic  moment  from  the  Magnus  moment,  although  this  is  not  required  for  the 
d  =  0  helical  motion  because  the  complex  angle  of  attack  is  also  zero.  Zero  spin  rate  is  obtained 
by  the  addition  of  a  second  angular  velocity  oriented  along  the  flight  body’s  longitudinal  axis 
that  is  equal  and  opposite  of  the  longitudinal  component  of  the  angular  velocity  required  to 
produce  the  rotation  about  the  axis  of  the  helix.  This  second  angular  velocity  has  no  effect  on 
the  transverse  angular  velocities  q  and  7. 


For  each  of  the  helical  motions,  the  transverse  aerodynamic  moment  in  the  nonrolling  frame  will 
be  periodic  in  time,  which  also  indicates  that  the  flow  field  will  be  periodic  in  time  when  viewed 
from  the  nonrolling  coordinate  frame.  However,  for  both  types  of  helical  motion,  a  steady  flow 
field  should  be  observed  when  examined  from  an  orthogonal  right-handed  coordinate  system  that 
has  its  x-axis  aligned  with  the  longitudinal  axis  of  the  body  and  its  z-axis  along  a  line  between 
the  body  CG  and  the  axis  of  rotation  of  the  helix.  The  transverse  aerodynamic  moments  in  the 
nonrolling  frame  can  be  transformed  to  the  coordinate  frame  where  the  steady  flow  field  exists, 
using  the  following  relations: 

q  =  Helical  Motion 

C„+iC.=(c„+icJe-“,  (15) 

d  =  0  Helical  Motion 

C„+iC.=(c„+iC.y““'.  (16) 

Here,  the  transverse  moments.  Cm  and  Cn,  refer  to  the  moments  about  the  y-  and  z-axes, 
respectively. 

Using  this  transformation,  the  transverse  moments  can  be  shown  to  have  the  following  form  for 
each  of  the  two  types  of  helical  motion. 


q  =  Helical  Motion 


C„  +iC„  =-C, 


Q£ 

vVy 


V 


Rq^ 

V 


(17) 
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d  =  0  Helical  Motion 


+iC„  =C. 


Q£ 

V 


R„Q 


V 


(18) 


The  resulting  expressions  for  the  transverse  moments  are  independent  of  time,  indieating  that  the 
flow  field  is  potentially  steady  as  well. 

Similar  expressions  for  the  individual  damping  foree  eoeffieients  ean  be  developed  using  the 
same  approaeh  as  applied  for  the  moment  eoeffieients. 


q  =  Helical  Motion 


Cy  +  iC^  -  iCfj 


vVy 


V 


Rq^ 

V 


«  =  0  Helical  Motion 


Cy  +iCz 


(19) 


(20) 


2,4  CG  Translation  Relations 


If  the  aerodynamie  eoeffieients  have  been  established  for  a  baseline  eonfiguration,  it  is  possible 
to  determine  the  aerodynamie  eoeffieients  for  the  identieal  eonfiguration  whieh  has  a  different 
axial  CG  loeation  using  the  CG  translation  relations  (12).  The  CG  translation  relations  for  the 
damping  eoeffieients  are  shown  in  equations  21-24.  In  the  eontext  of  the  eurrent  study,  the  CG 
translation  relations  ean  be  used  to  further  validate  the  results  by  eomparing  the  theoretieal 
variation  of  the  damping  eoeffieients  with  CG  loeation  with  the  CFD  predietions. 


C 


=  c 


(21) 


Cn.  -  Cn  +s^,gCN^  , 


+ScgCN^  , 


(22) 

(23) 


and 


C..  =C,^-s,g(C^^-C^J-s^gC 


(24) 


3.  Computational  Approach 


Computation  of  the  viseous  flow  field  about  the  flight  body  was  aecomplished  by  solving  the 
thin-layer  Navier-Stokes  equations  using  a  parabolized  Navier-Stokes  (PNS)  teehnique.  Beeause 
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the  computations  are  performed  in  a  noninertial  rotating  coordinate  frame,  the  governing 
equations  have  been  modified  to  include  the  body  force  terms  that  result  from  the  Coriolis  and 
centripedal  accelerations  in  the  rotating  coordinate  frame.  The  fluid  flow  relative  to  the  rotating 
coordinate  frame  does  not  vary  with  time,  allowing  the  steady  (nontime  varying)  Navier-Stokes 
equations  to  be  applied.  The  solution  of  the  steady  Navier-Stokes  equations  can  be  performed  at 
a  reasonable  computational  cost.  The  steady  thin-layer  Navier-Stokes  equations  are  shown  as 
follows: 


dE  5F  5G  .-v  1 

- 1 - 1 - 1"  H  +  H  —  — 

drj  d(^  Re 


0^  ^ 


(25) 


Here,  E ,  F,  and  G  are  the  inviscid  flux  vectors,  S  is  the  viscous  flux  vector,  and  S^,  are 

inviscid  and  viscous  source  terms  due  to  the  cylindrical  coordinate  formulation,  and  H  is  the 
source  term  containing  the  Coriolis  and  centrifugal  force  terms  that  result  from  the  rotating 
coordinate  frame.  Each  of  these  matrices  is  a  function  of  the  dependent  variables  represented  by 
the  vector  q^  =  {p,  pa,  pw,  pm,  e)  ,  where  p  and  e  are  the  density  and  the  total  energy  per  unit 
volume,  and  u,  v,  and  w,  are  the  velocity  components  in  axial,  circumferential,  and  normal 
directions.  The  flux  terms  are  shown  as  follows: 


where 


and 


pU 

pW 

/xrU  +  ^,,p 

F  =  - 
J 

paN  + 

pw\J 

PvV  +  7^p/r 

ycnvU 

yCWV  -1-  7j.P 

(e  +  p)U  J 

(e  +  p)V 

pW 

0 

puw  + 

H  =  i 
J 

pfx 

PvW  +  ^^p/r 

pwW  +  (^rP 

pfr 

(e-i-p)W 

/7u4  +/7vf^  +/7wf. 

U  =  u4 

V  =  / r  +  w;;^ 

W  =  u^,,  +v^^/r  +  w4',  ’ 


(26) 


(27) 
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(28) 


?7x 

rir=A-^^h) 

Cr  = 

J  =  l/[x^(^/i^r^  -^^r^)]. 

The  Coriolis  and  centrifugal  acceleration  terms,  due  to  the  rotating  coordinate  system,  which  are 
contained  in  the  source  term,  H ,  are  shown  as  follows: 

f  =  2Q X u  +  Qx  (q X  r).  (29) 

The  Coriolis  acceleration  is  a  function  of  the  angular  velocity  of  the  coordinate  frame  with 
respect  to  the  inertial  frame,  Q ,  and  the  fluid  velocity  vector,  u,  which  can  be  represented  by  the 
velocity  components,  u,  v,  and  w.  The  centripedal  acceleration  is  a  function  of  the  angular 
velocity  of  the  rotating  frame,  Q ,  and  the  displacement  vector,  R,  between  the  axis  of  rotation 
and  the  local  position  in  the  flow  field.  The  acceleration  vector,  f ,  can  be  written  in  terms  of  its 
components  along  the  x,  (f)  and  r  axes,  f^ ,  f^ ,  and  f^ . 


The  angular  velocities  for  the  three  different  motions  considered  in  this  study  are  written  in  the 
following.  The  angular  velocities  are  nondimensionalized  in  a  manner  consistent  with  the 
Navier-Stokes  equations. 


Looping  Motion 


d  =  0  Helical  Motion 


Q  = 


(  QD^ 

f  QD^ 

cos(pi^  - 

y^cc  J 

(  ^co  J 

sin  (j)  ij. 


Q  = 


V  ^00  y 


(30) 


(31) 


q  =  Helical  Motion 

Q  = 


( 

f  QD^ 

f  QD^ 

cos  R  + 

sin  cos 

V  ^00  ) 

1^00  ; 

1^00  J 

sin  P  sin  (j)  i^ . 


(32) 


The  pressure,  p,  can  be  related  to  the  dependent  variables  by  applying  the  ideal  gas  law: 


p  =  (/-!) 


e-— (u^  +  v^  +  w^ 
2  ^ 


(33) 
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The  turbulent  viscosity,  ju^,  which  appears  in  the  viscous  matrices,  was  computed  using  the 
Baldwin-Lomax  turbulence  model  {13). 

The  thin-layer  equations  are  solved  using  the  PNS  technique  of  Schiff  and  Steger  {14). 

Following  the  approach  of  Schiff  and  Steger,  the  governing  equations,  which  have  been  modified 
here  to  include  the  Coriolis  and  centrifugal  force  terms,  are  solved  using  a  conservative, 
approximately  factored,  implicit  finite-difference  numerical  algorithm  as  formulated  by  Beam 
and  Warming  {15).  Details  of  the  implementation  of  the  source  term  that  contains  the  Coriolis 
and  centrifugal  force  terms  are  given  in  reference  (5). 

The  computations  presented  here  were  performed  using  a  shock-fitting  procedure  {16).  This 
procedure  solves  the  five  Rankine-Hugoniot  jump  conditions,  two  geometric  shock  propagation 
conditions,  and  one  compatibility  equation  to  determine  the  values  of  the  five  dependent 
variables  immediately  behind  the  shock,  as  well  as  the  position  of  the  shock.  By  including  the 
implicit  part  of  the  source  term  due  to  the  rotating  coordinate  frame  in  the  circumferential 
inversion,  the  shock-fitting  procedure  can  be  used  without  modification,  as  long  as  the  free- 
stream  conditions  are  modified  to  account  for  the  rotating  coordinate  frame. 

At  the  body  surface,  no-slip,  constant  wall  temperature  boundary  conditions  were  applied.  For 
the  helical  motions,  the  tangential  velocity  at  the  body  surface  was  set  equal  to  the  local  velocity 
of  the  body  surface  due  to  the  solid-body  rotation  produced  by  an  angular  velocity  equal  and 
opposite  to  the  longitudinal  component  of  angular  velocity  shown  in  equations  31  and  32.  This 
is  required  to  produce  a  zero  spin  rate  p  in  the  nonrolling  coordinate  frame. 

An  initial  solution  for  the  PNS  marching  procedure  was  obtained  using  a  conical  step-back 
procedure  at  a  location  of  0.2  body  diameters  from  the  nosetip.  Although  the  perturbations  to  the 
flow  field  from  the  helical  motions  are  not  compatible  with  the  conical  flow  assumption,  the 
effect  on  the  solution  appears  to  be  small.  Moving  the  starting  plane  to  0.1  body  diameters  from 
the  nose  tip  resulted  in  <0.5%  variation  in  the  computed  damping  moments. 

The  computational  results  presented  here  were  obtained  using  a  grid  that  consisted  of  60  points 
between  the  body  and  the  shock.  In  the  circumferential  direction,  gridding  was  performed  over  a 
360°  sector  because  of  the  lack  of  symmetry  from  the  combination  of  angle  of  attack,  spin,  and 
coning  or  helical  motion.  Thirty-six  grid  points  were  used  in  the  circumferential  direction.  In 
longitudinal  direction,  78  marching  steps  were  utilized  for  each  body  diameter  of  length.  To 
ensure  adequate  grid  resolution  within  the  boundary  layer,  the  grid  spacing  at  the  body  surface 
was  adapted  to  maintain  nondimensional  boundary  layer  coordinate  y^  between  2  and  3  in 
accordance  with  previously  published  results  {1 7).  Grid  resolution  studies  for  the  9-caliber  body 
(middle  CG  location)  showed  <2%  variation  in  the  computed  pitch-damping  coefficient  sum 
when  the  grid  resolution  was  decreased  by  25%  in  each  of  the  three  coordinate  directions. 
Similarly,  when  the  grid  in  the  circumferential  and  marching  directions  was  doubled  and  the  grid 
in  the  radial  direction  was  increased  by  50%,  the  computed  pitch-damping  coefficient  sum  varied 
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by  only  1%.  Differences  in  the  damping  coefficient  on  the  coarse  and  fine  mesh,  relative  to 
the  baseline  grid,  were  <1%.  The  damping  coefficient  showed  a  larger  variability  on  a 

percentage  basis  on  the  coarse  and  fine  grid  relative  to  the  baseline  grid  (6%  and  3%, 
respectively),  although  absolute  differences  between  the  predicted  values  of  the  coefficient  were 
similar  to  the  differences  between  the  predicted  pitch-damping  coefficient  sum.  The  smaller 
magnitude  of  ,  relative  to  the  pitch-damping  coefficient  sum,  contributed  to  a  larger 

percentage  error.  The  computations  for  the  baseline  grid,  which  were  performed  using  a  Cray  Y- 
MP  supercomputer,  typically  required  <10  min  of  CPU  time  for  complete  calculation  over  a 
single  configuration. 


4.  Results 


Computations  of  the  individual  coefficients  that  comprise  the  pitch-damping  force  and  moment 
coefficient  sums  were  performed  using  helical  and  looping  motions.  The  computations  were 
performed  for  the  ANSR  series  of  bodies  shown  in  figure  1.  Results  for  length-to-diameter 
(L/D)  ratio  =  5,  7,  and  9  bodies  are  presented  here  for  flight  velocities  of  Mach  1.8  and  2.5. 

Three  different  CG  locations  were  considered  for  each  body  length,  as  shown  in  table  1 . 
Aerodynamic  range  tests  of  the  ANSR  (18)  were  used  to  benchmark  the  predicted  pitch-damping 
coefficient  sum. 


Table  1.  CG  location  for  ANSR  flight  bodies. 


Longitudinal  CG 
(Calibers  From  Nose) 

L/D 

Forward 

Middle 

Rearward 

5 

2.5 

3.0 

3.5 

7 

3.251 

4.037 

4.812 

9 

4.0 

5.038 

5.885 

The  effect  of  angular  rate  and  the  rotational  velocity  ratio  for  the  q  =  0  helical  motion  and  dr  =  0 
helical  motion  was  examined  for  the  L/D  =  9  ANSR  body  (forward  CG)  at  a  flight  velocity  of 
Mach  2.5.  The  force  and  moment  expansion  for  both  of  these  motions  indicates  a  linear 
variation  of  the  forces  and  moments  is  expected  with  angular  rate  and  rotational  velocity  ratio 
within  the  regime  where  linear  aerodynamic  theory  is  expected  to  be  valid.  The  predictions  were 
performed  for  the  longest  ANSR  body  because  the  longer  bodies  typically  exhibit  more 
nonlinear  aerodynamic  behavior.  Thus,  the  onset  of  nonlinear  behavior  of  the  aerodynamics 
with  angular  rate  or  angle  of  attack  would  be  more  evident,  and  the  limits  of  linear  aerodynamic 
behavior  could  be  more  readily  identified. 
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The  effect  of  angular  rate  for  the  q  =  0  motion  was  examined  by  performing  computations  at 
several  angular  velocities  but  at  a  fixed  total  angle  of  attack  of  2°.  This  means  that  as  the  angular 
velocity  was  increased,  the  radius  Ro  was  decreased  so  that  the  product  RqO  (and  the  total  angle 

of  attack)  was  constant.  Figure  8  shows  the  variation  of  the  aerodynamic  moment  Cm  with 
angular  velocity,  and  the  variation  is  seen  to  be  linear  across  the  range  of  angular  velocities  of 
interest  here. 


Figure  8.  Variation  of  aerodynamic  moment  Cm  with  angular  velocity,  q  =  0  helical  motion, 
M  =  2.5,  ANSR,  sin"'  S  =  2°. 


For  the  q  =  0  helical  motion,  the  sine  of  the  angle  of  attack,  S ,  varies  directly  with  the  rotational 
R  Q 

velocity  ratio,  ^  (see  equation  8).  Thus,  as  the  rotational  velocity  ratio  increases  toward  1, 

the  angle  of  attack  of  the  vehicle  also  increases  for  the  q  =  0  helical  motion.  (In  the  limit  as 
R  Q 

— =  1 ,  for  q  =  0  helical  motion,  the  CG  performs  a  circular  looping  motion  with  the  body  at 

90°  angle  of  attack.)  The  desire  here  is  to  define  the  limits  where  the  forces  and  moments  are 

R  Q 

expected  to  vary  linearly  with  S  or  — .  Computations  of  q  =  0  helical  motion  were 

QD 

performed  with  the  angular  rate  held  fixed  ( =  0.01 ),  while  the  helix  radius  was  varied  to 
produce  the  desired  angle  of  attack.  The  aerodynamic  moment  Cm,  shown  in  figure  9,  shows  a 
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Figure  9.  Variation  of  aerodynamic  moment  Cm  with  sine  of  the  angle  of  attack,  q  =  0  helical 


motion,  M  =  2.5,  ANSR, 


QD 

V 


linear  variation  with  S  up  to  angles  of  attack  of  ~7°.  This  behavior  is  similar  to  that  observed 
for  coning  motion  (5). 

For  the  dr  =  0  helical  motion,  the  effect  of  angular  rate  was  examined  for  rotational  velocity 
QR 

ratio,  =  sin  =  0.0349  .  Again,  as  the  angular  velocity  was  increased,  the  radius  Ro  was 

decreased  so  that  the  product  RqO  (and  sin  )  was  constant.  The  variation  of  the  aerodynamic 
moment  with  angular  velocity,  shown  in  figure  10,  is  seen  to  be  linear  across  the  range  of 
angular  velocities  of  interest  here. 


The  effect  of  the  rotational  velocity  ratio  ^  °  =  sin  on  the  forces  and  moments  for  the  dr  =  0 

helical  motion  was  also  examined,  and  results  are  shown  in  figure  1 1 .  These  predictions  were 

QD 

made  for  a  constant  angular  velocity,  ( =  0.01 ).  The  aerodynamic  moment  Cm  is  seen  to  be 

linear  with  rotational  velocity  ratios  from  0  to  1 .  At  first  glance,  a  nonlinear  variation  might  be 
expected  since  nonlinear  variations  with  S  were  seen  for  the  q  =  0  helical  motion.  It  is  noted 
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Figure  10.  Variation  of  aerodynamic  moment  with  angular  velocity,  d  -  0  helical  motion, 
M  =  2.5,  ANSR,  p  =  2°. 
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QR 

that  the  parameter  ^  °  =  sin  P  is  not  assoeiated  with  the  angle  of  attaek  for  this  motion 

beeause  the  dr  =  0  helical  motion  is  performed  at  0°  angle  of  attack.  The  parameter,  sin  P,  is 
simply  the  rotational  velocity  ratio  which  for  the  dr  =  0  helical  motion  reflects  the  aspect  ratio  of 
the  helix.  The  limiting  case  of  a  rotational  velocity  ratio  of  1  corresponds  to  a  circular  looping 
motion  discussed  previously. 

Using  these  results  as  a  guide,  predictions  of  the  pitch-damping  coefficients  were  made  for  each 
of  the  three  body  lengths  and  three  CGs  at  Mach  1 .8  and  2.5.  Predictions  were  made  at  an 
QD  QR 

angular  velocity  of  =  0.01  and  a  rotational  velocity  ratio,  ^  °  =  0.0349,  where  the  linear 
aerodynamic  theory  has  been  shown  to  be  valid. 

Figures  12-14  show  predictions  of  the  pitch-damping  moment  coefficients,  and  ,  as  a 

function  of  the  CG  position  for  the  L/D  =  5,  7,  and  9  bodies  at  Mach  2.5.  Predictions  obtained 
using  the  PNS  computational  approach  are  shown,  along  with  results  obtained  with  slender  body 

theory.  Both  PNS  and  slender  body  predictions  show  that  the  pitch-damping  coefficient  is 

larger  than  for  all  the  CG  positions  examined.  The  PNS  prediction  of  the  damping 

coefficient,  ,  is  nearly  zero  at  the  rearward  CG  for  each  of  the  body  lengths.  The  trends 

shown  by  the  slender  body  results  are  similar  to  the  PNS  results,  although  the  slender  body 
predictions  are  generally  lower  in  magnitude.  The  PNS  predictions  of  the  variation  of  Cj^^  with 

CG  position  shown  in  figures  12-14  agree  with  the  theoretical  variation  (equation  24)  to  within 
0.25%  across  the  range  of  CG  positions  examined.  The  PNS  predictions  of  the  variation  of  Cj^^ 

are  essentially  identical  to  the  theoretical  variation  (equation  23)  because  the  governing 
equations,  boundary  conditions,  and  flow  field  are  independent  of  the  axial  CG  location. 

Figures  15-17  show  similar  predictions  of  the  pitch-damping  force  coefficients,  Cj.^^  and  C^^  , 

as  a  function  of  CG  position  for  the  L/D  =  5,  7,  and  9  bodies  at  Mach  2.5.  The  PNS  predictions 
of  the  pitch-damping  force  coefficients  show  a  similar  trend  to  the  predictions  made  with  slender 
body  theory,  though  the  slender  body  theory  results  are  larger  in  magnitude.  The  agreement  does 
seem  to  improve  slightly  as  the  body  length  increases. 

Because  ballistic  bodies  traverse  a  nearly  rectilinear  path,  it  is  impossible  to  extract  the 
individual  components  of  the  damping  coefficients  from  the  experimental  data.  Thus,  no 
comparison  of  the  predicted  individual  coefficients  could  be  made  with  experimental  data. 
However,  the  individual  coefficients  can  be  summed  to  produce  the  pitch-damping  coefficient 
sum,  and  the  results  compared  with  the  coning  motion  results  and  with  experiment.  Figures 
1 8-20  show  comparisons  of  the  pitch-damping  coefficient  sum  (obtained  by  summing  the 
individual  coefficients)  with  experimental  data  and  slender  body  theory  results.  The  results 
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FICIENT  f  DAMPING  COEFFICIENT 


:e  12.  Variation  of  damping  moment  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  5. 


CG  -  CALIBERS  FROM  NOSE 


Figure  13.  Variation  of  damping  moment  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  7. 
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DAMPING  FORCE  COEFFICIENT  f  DAMPING  COEFFICIENT 


-100.0 


MACH  NUMBER  =  2.5 
L/D=9. 


CG  -  CALIBERS  FROM  NOSE 


re  14.  Variation  of  damping  moment  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  9. 


CG  -  CALIBERS  FROM  NOSE 


Figure  15.  Variation  of  damping  force  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  5. 
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DAMPING  FORCE  COEFFICIENT  f  DAMPING  FORCE  COEFFICIENT 


16.  Variation  of  damping  force  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  7. 


CG  -  CALIBERS  FROM  NOSE 


Figure  17.  Variation  of  damping  force  coefficients  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  9. 
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Figure  18.  Variation  of  damping  moment  coefficient  sum  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  5. 


CG  -  CALIBERS  FROM  NOSE 


Figure  19.  Variation  of  damping  moment  coefficient  sum  with  CG  location,  M  =  2.5,  ANSR, 
L/D  =  7. 


23 


CG  -  CALIBERS  FROM  NOSE 


Figure  20.  Variation  of  damping  moment  coefficient  sum  with  CG  location,  M  =  2.5,  ANSR, 

L/D  =  9. 

obtained  by  summing  the  predicted  coefficients,  and  ,  are  within  1.0%  or  less  of  the 
pitch-damping  moment  coefficient  sum  predicted  using  coning  motion. 

To  provide  additional  validation  of  the  current  approach,  the  predictions  were  also  compared 
with  the  nonlinear  unsteady  potential  equation  results  obtained  by  Devan  {19)  for  the  L/D  =  5 
ANSR  configuration.  The  variation  in  the  two  damping  coefficients  with  Mach  number  is  shown 
in  figure  21.  In  general,  the  results  show  modest  variations  in  the  predicted  coefficients  across 
the  Mach  number  range  except  for  the  slender  body  predictions  which  are  independent  of  Mach 
number.  The  nonlinear  potential  solution  of  Devan  shows  good  agreement  with  the  PNS 
predictions  with  a  similar  variation  with  Mach  number  predicted.  (It  is  noted  that  Devan  did  not 
explicitly  determine  the  damping  coefficient  ,  but  instead  computed  and 

(Cjjj  +  Cjjj .  )  from  which  Cj^  can  be  obtained.)  The  approach  of  Devan  overpredicts  the  two 

CC  Q 

damping  coefficients  by  7%  for  and  ~25%  for  ,  compared  with  the  PNS  results.  The 
predicted  value  of  from  slender  body  theory  differs  by  as  much  as  25%  with  the  PNS 

predictions.  The  slender  body  theory  predicts  a  nearly  zero  value  for  the  damping  coefficient 
,  while  the  PNS  predictions  show  this  coefficient  is  ~  25%-40%  of  the  value  of  at  the 

middle  CG  location. 
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Figure  21.  Variation  of  the  pitch-damping  moment  coefficients  with  Mach  number,  ANSR, 

L/D  =  5,  middle  CG  location. 

A  comparison  of  the  Mach  number  variation  of  the  pitch-damping  coefficient  sum  for  the  various 
approaches  with  aerodynamic  range  data  is  shown  in  figure  22.  Compared  with  the  slender  body 
theory  and  the  nonlinear  unsteady  potential  equation  results  of  Devan,  the  PNS  results  show  the 
best  agreement  with  the  experimental  data.  The  approach  of  Devan  shows  generally  a  similar 
variation  with  Mach  number  as  the  PNS  approach  and  the  experimental  data,  but  appears  to  show 
a  slight  overprediction  (-'10%)  across  the  Mach  number  range.  This  overprediction  is  consistent 
with  the  overprediction  seen  previously  for  the  individual  coefficients.  The  slender  body  theory 
results,  which  exhibit  no  dependence  on  Mach  number,  consistently  underpredicts  the  pitch¬ 
damping  sum  by  >30%. 

The  distribution  of  the  pitch-damping  force  and  moment  coefficients  over  the  ANSR  L/D  =  5 
body  is  shown  in  figures  23  and  24.  The  force  coefficient  is  positive  along  the  body  length 

with  most  of  the  force  being  generated  at  the  aft  end  of  the  body.  As  a  result,  the  moment 
coefficient  also  shows  its  largest  contributions  from  the  aft  portion  of  the  body.  On  the 

other  hand,  for  the  force  coefficient,  ,  both  the  nose  and  cylindrical  afterbody  produce 

contributions  to  the  force  coefficient  which  are  similar  in  magnitude  but  opposite  in  sign.  This 
essentially  imposes  a  couple  on  the  body  whose  magnitude  is  reflected  in  the  moment 
coefficient,  Cjjj  . 
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Figure  22.  Variation  of  the  piteh-damping  moment  eoeffieient  sum  with  Maeh  number, 
ANSR,  L/D  =  5,  middle  CG  loeation. 


Figure  23.  Development  of  piteh-damping  foree  eoeffieients  over  ANSR  body,  M  =  2.5, 
L/D  =  5,  middle  CG  position. 
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Figure  24.  Development  of  pitch-damping  moment  coefficients  over  ANSR  body,  M  =  2.5, 
L/D  =  5,  middle  CG  position. 


5.  Conclusion 


A  method  for  determining  the  individual  force  and  moment  coefficients  that  comprise  the  pitch¬ 
damping  force  and  moment  coefficient  sum  has  been  presented,  along  with  sample  results  for  a 
family  of  axisymmetric  projectile  geometries.  The  sum  of  the  two  individual  coefficients  is 
identical  to  previous  predictions  of  the  pitch-damping  coefficient  sums  obtained  using  coning 
motion  and  is  in  excellent  agreement  with  experimental  data.  Additional  validation  of  the 
approach  is  obtained  through  comparisons  with  prior  results  from  numerical  solutions  of  the 
nonlinear  unsteady  potential  equation.  The  individual  coefficients  show  qualitative  agreement 
with  results  obtained  using  the  more  approximate  slender  body  theory.  The  method  provides  an 
efficient  means  of  determining  the  individual  coefficients  which  may  be  quite  difficult  to  obtain 
using  experimental  means. 
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List  of  Abbreviations  and  Symbols 


Freestream  speed  of  sound 


Pitehing  moment  eoeffieient, 


Pitehing  moment  eoeffieient  slope  with  respect  to  angle  of  attack, 


Pitch-damping  moment  coefficient  slope, 


Pitch-damping  moment  coefficient  slope. 


Cmq  +  Pitch-damping  moment  coefficient 

C„  Side  moment  coefficient 


Magnus  moment  coefficient. 


d{  ^]da 


Normal  force  coefficient. 


Normal  force  coefficient  slope  with  respect  to  angle  of  attack, - ^ 

da 


Pitch-damping  force  coefficient  slope. 


Pitch-damping  force  slope. 


Pitch-damping  force  coefficient 
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D 


Reference  diameter 


tF,G 

f  f  f 

-^X  ’  ’  -^r 

F 

H 

He 

J 

£ 

L 

M 

M 

P 

P 

q,r 

r 

R 

Ro 

Re 

s 

^cg 

S.f 

s 

Se 

t 


Total  energy  per  unit  volume 

Flux  vectors  in  transformed  coordinates 

Component  of  acceleration  vector  f  in  the  x ,  ^ ,  and  r  directions 
Force 

Source  term  inNavier-Stokes  equations  due  rotating  coordinate  frame 
Inviscid  source  term  due  to  cylindrical  coordinate  formulation 
Jacobian 

Characteristic  length,  typically  D 

Body  length 

Moment 

Freestream  Mach  number 

Pressure,  as  used  in  Navier-Stokes  equations 

Spin  rate,  as  used  in  roll  equations 

Transverse  angular  rates  of  body  along  y  and  z  directions 

Radial  coordinate  as  used  in  Navier-Stokes  equations 

Displacement  vector  between  the  axis  of  rotation  and  location  in  flow  field 

Helix  radius  or  radius  of  circular  arc  for  looping  motion 

Reynolds  number,  a^p^D/  ju^ 

Distance  downrange 

Center  of  gravity  shift  in  calibers 


Reference  area. 


Viscous  flux  vector 


Additional  viscous  terms  due  to  cylindrical  coordinate  formulation 


Time 


32 


u,  V,  w  Velocity  components  in  the  three  directions 

U  Vehicle  velocity  alone  helix  axis 

Y  Magnitude  of  freestream  velocity 

V^  Freestream  velocity  vector 

X,  y,z  Axial,  horizontal,  and  vertical  coordinates 

X g ,  y g ,  Zg  Three-dimensional  coordinate  in  earth- fixed  system 

Xgg  Axial  location  of  center  of  gravity 

y^  Nondimensional  boundary  layer  coordinate 

Greek  Symbols 

a  Angle  of  attack 

d  Angular  rate  associated  with  angle  of  attack 

[5  Angle  between  helix  axis  and  longitudinal  body  axis  for  d  =  0  helical  motion 

Y  Ratio  of  specific  heats,  in  Navier-Stokes  equations 

Y  Cosine  of  total  angle  of  attack,  as  used  in  the  aerodynamic  force  and  moment 
equations 

5  Sine  of  the  total  angle  of  attack 

0  Angular  orientation  of  body  relative  to  inertial  frame  for  looping  motion 

//,//(  Laminar  and  turbulent  viscosity 

Jl  Complex  transverse  angular  velocity 

7,  Transformed  coordinates  in  the  Navier-Stokes  equations 

^  Complex  angle  of  attack 

p  Density 

(j)  Circumferential  coordinate,  measured  from  the  vertical  axis 

Q  Angular  rate  associated  with  looping  and  helical  motions 
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Superscripts 

C) 

o' 

(') 

Subscripts 


Rate  of  change  with  respect  to  time 

Rate  of  change  with  respect  to  space 
Referenced  to  nonrolling  coordinate  frame 


00 


Quantity  evaluated  at  freestream  conditions 
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